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We investigate the relevant gluon-momentum region for confinement in lattice QCD on 16 4 at 
/3=5.7, 5.8 and 6.0, based on the Fourier expansion. We find that the string tension a, i.e., the 
confining force, is almost unchanged even after removing the high-momentum gluon component 
(N. above 1.5GeV in the Landau gauge. In fact, the confinement property originates from the low- 

momentum gluon component below 1.5GeV, which is the upper limit to contribute to a. In the 
relevant region, smaller gluon momentum component is more important for confinement. Next, we 
develop a manifestly gauge-covariant expansion of the QCD operator such as the Wilson loop, using 
5-H ' the eigen-mode of the QCD Dirac operator If) = 7 M D M . With this method, we perform a direct 

analysis of the correlation between confinement and chiral symmetry breaking in lattice QCD on 
6 4 at /3=5.6. As a remarkable fact, the confinement force is almost unchanged even after removing 
the low-lying Dirac modes, which are responsible to chiral symmetry breaking. This indicates that 
one-to-one correspondence does not hold for between confinement and chiral symmetry breaking in 
QCD. 
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I. INTRODUCTION 

ID; 

Nowadays, quantum chromodynamics (QCD) has been established as the fundamental gauge theory of the strong 
interaction. However, nonperturbative properties of low-energy QCD such as color confinement and chiral symmetry 
breaking [l[ are not yet well understood, which gives one of the most difficult problems in theoretical physics. The 
nonperturbative QCD has been studied in lattice QCD [2|44[ and various analytical frameworks [Bl-flll]. 

In this paper, using lattice QCD, we research for the origin of color confinement in terms of the relevant gluon- 
momentum component based on the Fourier expansion in the Landau gauge [l2j . We also investigate the correspon- 
dence between color confinement and chiral symmetry breaking using the Dirac-mode expansion in a gauge-invariant 
manner. 



II. RELEVANT REGION OF GLUONIC MOMENTUM FOR COLOR CONFINEMENT 

Many theoretical physicists consider that confinement phenomenon is brought by low-energy region of QCD, because 
of the strong QCD coupling in the infrared (IR) region. However, at the quantitative level, it is difficult to state the 
"relevant energy region" for confinement directly from QCD. Since nonperturbative phenomena are mainly brought by 
gluon dynamics, the key question here is "what is the relevant gluon-momentum region responsible for confinement?" 

In this section, to get the answer, we study quantitative lattice-QCD analysis for the relevant gluon-momentum 
region for color confinement [l2[, based on the Fourier expansion of the link- variable. Our method consists of the 
following five steps. 

Step 1. Generation of link-variable in the Landau gauge 
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We generate a gauge configuration on a L 4 lattice with the lattice spacing a by the lattice-QCD Monte Carlo 
method under space-time periodic boundary conditions. Here, we consider the link- variable U^{x) = e %a 9 A A x ) fixed 
in the Landau gauge, where the fluctuation from gauge degrees of freedom is strongly suppressed, owing to the global 
suppression of gluon- field fluctuations [f3j . 

Step 2. Four- dimensional discrete Fourier transformation 

By the discrete Fourier transformation, we define the "momentum-space link- variable" , 

U^(p) = ^r^— V U^x) exp^vPvXv), (f ) 

jVsito 

with the lattice-site number N s n e . The momentum-space lattice spacing is given by a p = 2ir/(La). 
Step 3. Imposing a cut in the momentum space 

We impose a cut on U fJl (p) in the momentum space, as shown in Fig.f . Outside the cut, we replace U^p) by the free- 
field link-variable, U { J ee (p) = jf— Yl x ex P(*S ! y? 3 ^ ;l; ^) = ^po, corresponding to U^x) = I. Then, the momentum-space 
link- variable U^{p) with the cut is defined as 

T~T K (n\ = / ^niv) (inside cut) 

" W \ U* cc (p) = S p0 (outside cut). [ ) 




FIG. 1: A schematic figure of the UV cut Auv and the IR cut Air on momentum-space lattice, with the lattice spacing 
a p = 2ir/(La). The momentum-space link- variable L^t(p) is replaced by the free-field link- variable U^ cc (p) = S p o in the shaded 
cut regions. 



Step 4- Inverse Fourier transformation 

To return to the coordinate space, we carry out the inverse Fourier transformation as 

W = E ^» expHE.^)- (3) 
v 

Since this U'^x) is not an SU(3) matrix, we project it onto an SU(3) element U^{x) by maximizing ReTr[J7^(x)'l'{7^(a;)]. 
Such a projection is often used in lattice QCD algorithms. By this projection, we obtain the coordinate-space link- 
variable U^(x) with the cut, which is an SU(3) matrix and has the maximal overlap to U'^x). 

Step 5. Calculation of physical quantities 

Using the cut link- variable U^(x), instead of Ufj,(x), we calculate physical quantities as the expectation value in the 
same way as original lattice QCD. 



3 



With this method in lattice-QCD framework, we quantitatively determine the relevant energy scale of color con- 
finement, through the analyses of the QQ potential. The lattice QCD Monte Carlo simulations are performed on 16 4 
lattice at /3=5.7, 5.8 and 6.0 at the quenched level (T^ |. 
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FIG. 2: (a) The QQ potential V(R) with the IR cut Air plotted against the inter-quark distance R. (b) The QQ potential with 
the UV cut Auv- Lattice QCD calculations are done on 16 4 lattice with j3 = 6.0, i.e., a ~ O.lOfm and a p = 2ir/(La) ~ 0.77GeV 
The dashed line is the original QQ potential in lattice QCD. 



Figure 2 (a) and (b) show the QQ potential V(R) with the IR cutoff Air, and the UV cutoff Auv, respectively. We 
get the following lattice-QCD results on the role of gluon momentum components. 

• By the IR cutoff Air, as shown in Fig. 2(a), the Coulomb potential seems to be unchanged, but the confinement 
potential is largely reduced fl2l ]. 

• By the UV cutoff Auv, a s shown in Fig. 2(b), the Coulomb potential is largely reduced, but the confinement 
potential is almost unchanged [12| . 
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FIG. 3: (a) The AiR-dependence and (b) the Auv-dependence of the string tension a. The string tension is obtained from 
the asymptotic slope of the QQ potential V(R) with the IR/UV cutoff in lattice QCD on 16 4 at j3 =5.7, 5.8 and 6.0 [3. The 
horizontal error-bar denotes the range from the discrete momentum, while the vertical error-bar is the statistical one. The 
broken line denotes the original value of a ~ 0.89GeV/fm. 



Figure 3 shows the AiR,/Auv-dcpcndcncc of the string tension a obtained from the asymptotic slope of the QQ 
potential V(R) with the IR/UV cutoff. Note that the ordinary QCD system without the cutoff corresponds to Air = 
and Auv = +oo. As shown in Fig. 3(a), the string tension is significantly reduced by the IR-cutoff Air even for small 
values of Air,, and smaller gluon-momentum component seems to be more important for confinement. 

As a remarkable fact, the string tension is almost unchanged even after cutting off the high-momentum gluon 
component above I.5GcV, as shown in Fig. 3(b) fl2j | - In fact, the confinement property originates from the low- 
momentum gluon component below 1.5GcV, which is the upper limit to contribute to a. Note here that the relevant 



4 



region \p\ < 1.5GcV for the confinement is only a small part of the total four-dimensional Brillouin zone of the gluon 
field, — ir/a < p^ < ir/a (yU=l,2,3,4). For example, at /3 = 6.0 (i.e., a ~ O.lfm), the relevant region \p\ < 1.5GeV is 
less than 0.2% in the total Brillouin zone, as shown in Fig. 4. 
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FIG. 4: The relevant gluon-momentum region \p\ < 1.5GeV for confinement in the Brillouin zone (BZ) —n/a < p M < 
7r/a(~ 6GeV) at /3 = 6.0. The relevant region is less than 0.2% in the four-dimensional BZ. In the relevant region, smaller 
gluon-momentum component is more important for confinement. 

With the same method, we find also the relevant role of low-momentum gluons to chiral symmetry breaking in 
lattice QCD [Til ]. However, the response pattern against the cut is rather different from that on confinement, and 
higher-momentum gluons also contribute to the chiral condensate. 

III. GAUGE-INVARIANT FORMALISM WITH DIRAC-MODE EXPANSION: 
A Direct Investigation of Correlation between Confinement and Chiral Symmetry Breaking 

Next, we newly develop a manifestly gauge-covariant expansion of the QCD operator such as the Wilson loop, using 
the eigen-mode of the QCD Dirac operator p = j^D^, and investigate the relation between confinement and chiral 
symmetry breaking. 

A. Gauge Covariant Expansion in QCD instead of Fourier Expansion 

The previous method is based on the Fourier expansion, i.e., the eigen-mode expansion of the momentum operator 
p^. Because of the commutable nature of [p^Sp"] = 0, all the momentum p^ can be simultaneously diagonalized. 
which is one of the strong merits of the Fourier expansion. Also it keeps Lorcntz covariance. 

However, the Fourier expansion does not keep gauge invariance in gauge theories. Therefore, for the use of the 
Fourier expansion in QCD, one has to select a suitable gauge such as the Landau gauge, where the gauge-field 
fluctuation is strongly suppressed in Euclidean QCD. 

As a next challenge, we consider a gauge-invariant method, using a gauge-covariant expansion in QCD instead of 
the Fourier expansion. In fact, we consider a generalization of the Fourier expansion or an alternative expansion with 
keeping the gauge symmetry. 

A straight generalization is to use the covariant derivative operator instead of the derivative operator d^. 
However, due to the non-commutable nature of [D^,D U ] ^ 0, one cannot diagonalize all the covariant derivative £) M 
(/j = 1,2,3,4) simultaneously, but only one of them can be diagonalized. For example, the eigen-mode expansion 
of -D4 keeps gauge covariance and is rather interesting, but this type of the expansion inevitably breaks the Lorcntz 
covariance. Then, we consider the eigen-mode expansion of the Dirac operator p = ^^D^ or D 2 = D^D^ : since such 
an expansion keeps both gauge symmetry and Lorentz covariance. 

In particular, the Dirac-mode expansion is rather interesting because the Dirac operator p directly connects with 
chiral symmetry breaking via the Banks-Casher relation Q and its zero modes are related to the topological charge 
via the Atiyah-Singer index theorem Here, we mainly consider the manifestly gauge-invariant new method using 
the Dirac-mode expansion. Thus, the Dirac-mode expansion has some important merits. 

• The Dirac-mode expansion method manifestly keeps both gauge and Lorentz invariance. 

• Each QCD phenomenon can be directly investigated in terms of chiral symmetry breaking. 
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B. Confinement and Chiral Symmetry Breaking: 
Can we expect One-to-one Correspondence between them in QCD? 



In particular, it is rather interesting and important to examine the correlation between confinement and chiral 
symmetry breaking 0, H, [HI HH, [lj|, since the direct relation is not yet shown between them in QCD. The strong 
correlation between them has been suggested by the simultaneous phase transitions of deconfinement and chiral 
restoration in lattice QCD both at finite temperature [3] and in a small- volume box 

The close relation between confinement and chiral symmetry breaking has been also suggested in terms of the 
monopole/vortex degrees of freedom [l5|, which topologically appears in QCD by taking the maximally 
Abelian/center gauge [El. Il7l. [l8j . For example, by removing the monopoles, confinement and chiral symmetry breaking 
are simultaneously lost in lattice QCD [15j . as schematically shown in Fig. 5. This indicates an important role of the 
monopolc to both confinement and chiral symmetry breaking, and these two nonperturbative QCD phenomena seem 
to be related via the monopole. 
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FIG. 5: An illustration of the relevant role of monopoles to nonperturbative QCD. In the maximally Abelian gauge, QCD 
becomes Abelian- like due to the large off-diagonal gluon mass of about lGeV [13] , and there appears a global network of the 
monopole current [l7l . Il8j . By the Hodge decomposition, the QCD system can be divided into the monopole part and the 
photon part. The monopole part has confinement [ll|, chiral symmetry breaking [llj and instantons [2(J, while the photon 
part does not have all of them. 



However, as a possibility, removing the monopoles may be "too fatal" for most nonperturbative properties. If this 
is the case, nonperturbative QCD phenomena are simultaneously lost by their cut. 

In fact, if only the relevant ingredient of chiral symmetry breaking is carefully removed, how will be confinement? 
To get the answer, we perform a direct investigation between confinement and chiral symmetry breaking, using the 
Dirac-mode expansion. 
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C. Eigen-mode of Dirac Operator in Lattice QCD 

In lattice QCD with spacing a, the Dirac operator p = ■y^D^ is expressed with U^{x) as 

I 4 

P x >y = Ya ^ 7m PvW^+hv ~ U -»( x ) S x-!i,y\ < ( 4 ) 

where U-^x) = UUx — jX). In the use of hermite 7- matrix 7^ = 7^, _p is anti-hermitc and satisfies ^ x = —p XtV . 
The normalized eigen-state |n) of the Dirac operator p is introduced as 

p\n) = i\ n \n) (5) 

with A„ 6 R. Because of {75, p} = 0, the state 75 |n) is also an eigen-state of p with the eigenvalue —i\ n - The Dirac 
cigenfunction ip n (x) = (x\n) obeys pip n (x) = i\ n i() n (x), and its explicit form of the eigenvalue equation in lattice 
QCD is 



I 4 

2~ X! 7/4^ 0*0 fa + A) ^ U-^x^nix - ji)] = i\ n ip n (x). (6) 

The Dirac eigenfunction ip n (x) can be numerically obtained in lattice QCD, besides a phase factor. 
According to U^(x) —> V(x)U IJl (x)V^ (x + fi), the gauge transformation of ip n (x) is found to be 

^ n {x)^V{x)ip n {x), (7) 

which is the same as that of the quark field. To be strict, for the Dirac eigenfunction, there can appear an irrelevant 
n-dependent global phase factor as e liPn ^ , according to the arbitrariness of the definition of ip n (x) . 

Note that the quark condensate (qq), the order parameter of chiral symmetry breaking, is given by the zero- 
eigenvalue density p(0) of the Dirac operator, via the Banks-Casher relation Q, 

{qq) = - Um Jim 7^(0). (8) 

m— >U V — foo 

Here, p(X) = y ^2 n ($(^ ~ ^n)) i s the spectral density of the Dirac operator. Also, the zero-mode number asymmetry 

of the Dirac operator p is equal to the topological charge (the instanton number) Q = J d 4 x Tr (G^G^), which 
is known as the Atiyah-Singer index theorem, Index(p)=Q Q. 



D. Operator Formalism in Lattice QCD 

To keep the gauge symmetry, careful treatments are necessary, since naive approximations may break the gauge 
symmetry. Here, we take the "operator formalism" , as explained below. 
We define the link-variable operator [/„ by the matrix element of 

(x\Up\y) = U ll {x)5 x+ p, >y . (9) 

The Wilson-loop operator W is defined as the product of XJ^ along a rectangular loop, 

N 

w = H% k = u l , 1 %2---% N - (10) 



fc=l 



For arbitrary loops, one finds Ylk=i Afe = 0- We note that the functional trace of the Wilson-loop operator W is 
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proportional to the ordinary vacuum expectation value (W) of the Wilson loop: 
Tr W ~ tTj2(x\W\x) =trJ2(x\U^U fl 2---U m \x) 

X X 

= tr Y (xi\U^ 1 \x 2 )(x 2 \U fl2 \x s ,)(x 3 \U l _ l3 \x 4 ) ■ ■ ■ (xTvl^jvki) 

Xx,X2,--- ,X N 

1 N-l 

= tT ^2( x \U^\x +fn)(x + Ail^Jx + ^2fi k ) ■ ■ ■ (x + Y MU m \x) 

x k=l k=l 

2 N-l 

= ^^{U^{x)U^{x + jli)U^{x + ^2(i k )---U llN {x+ Y Afe)} 

a- k=l k=l 

= (W)-Tr 1. (11) 

Here, "Tr" denotes the functional trace, and "tr" the trace over SU(3) color index. 

The Dirac-mode matrix element of the link- variable operator can be expressed with tp n (x) : 

(m\U\n) = 'Y^(m\x)(x\U fl \x + jl)(x + fi\n) = Y^L^W^Hnix + /t). (12) 

X X 

Although the total number of the matrix element is very huge, the matrix element is calculable and gauge invariant, 
apart from an irrelevant phase factor. Using the gauge transformation ([7]), we find the gauge transformation of the 
matrix element as 

(m\Ufj,\n) = ^2^ m (x)U^(x)^ n (x + fi) 

X 

-> J2 ^m^)V\x) ■ V{X)U^X)V\X + £) ■ V(x + + /i) 

X 

= Y,^l(x)Mx)iPn(x + ti = HU tl \n). (13) 

x 

To be strict, there appears an n-dependent global phase factor, corresponding to the arbitrariness of the phase in the 
basis \n). However, this phase factor cancels as e~ vpn e VPn = 1 between \n) and (n\, and does not appear for QCD 
physical quantities including the Wilson loop. 



E. Dirac-mode Expansion and Projection 

From the completeness of the Dirac-mode basis, J2 n \ n )( n \ = 1' we S et ^ = EmEn \m) (m\0\n) {n\ for arbitrary 
operators. Based on this relation, the Dirac-mode expansion and projection can be defined. We define the projection 
operator P which restricts the Dirac-mode space, 

P=J2\n){n\, (14) 

where A denotes arbitrary set of Dirac modes. In P, the arbitrary phase cancels between \n) and (n\. One finds 
P 2 = P and P^ = P. The typical projections are IR-cut and UV-cut of the Dirac modes: 

p m = y i n >H Pw = Y i n >H ( 15 ) 

A„|>Am |A„|<Auv 

Using the projection operator P, we define the Dirac-mode projected link-variable operator, 

U£ ee P%P =YJ2 \m){m\%\n){n\. (16) 

meA neA 

During this projection, there appears some nonlocality in general, but it would not be important for the argument 
of large-distance properties such as confinement. From the Wilson-loop operator W = IIaLi ' wc define the 
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Dirac-modc projected Wilson-loop operator, 

N 

w p = n k = k k ■ ■ ■ C = PK pk p- p u» N p 



fc=l 

2J |ni)(ni|{7 (Ul |n.2)(n.2|J7 M2 |n 3 ) • • • (njv|J7 MJV |njv+i)(njv+i|. (17) 



Then, we obtain the functional trace of the Dirac-mode projected Wilson-loop operator, 

N 

Tr W p = Tr J] tf£ = Tr • • • ^ P v = Tr PU^PU^P • • ■ Pl^P 

fe=i 

= tr ^ {n\\U^\n2){n2\U^ 2 \ns) ■■ ■ {n N \U m \ni), (18) 

which is manifestly gauge invariant. Here, the arbitrary phase factor cancels between \rik) and (nk\- Its gauge 
invariance is also numerically checked in the lattice QCD Monte Carlo calculation. 

From Tr W p (R,T) corresponding to the RxT rectangular loop, we define the Dirac-mode projected inter-quark 
potential, 

V P (R) ee - lim iln{Tr W P (R, T)}, (19) 

T— too 1 

which is also manifestly gauge-invariant. In the unprojected case of P = 1, the ordinary inter-quark potential is 
obtained apart from an irrelevant constant, 

V(R) = - lim -ln{TrW{R,T)} = - lim —\n(W(R,T)} + irrelevant const., (20) 
because of Tr W — (W) ■ Tr 1, as was derived in Eq. fTTl) . 

F. Analysis of Confinement in terms of Dirac Modes in QCD 

We consider various projection space A in the Dirac-mode space, e.g., IR-cut or UV-cut of Dirac modes. With 
this Dirac-mode expansion and projection formalism, we calculate the Dirac-mode projected inter-quark potential 
V P (R) in a gauge- invariant manner. In particular, using IR-cut of the Dirac modes, we directly investigate the role 
of low- lying Dirac modes to confinement, since the low- lying modes are responsible to chiral symmetry breaking. 

As a technical difficulty of this formalism, we have to deal with huge dimensional matrices and their products. 
Actually, the total matrix dimension of (m\U^\n) is (Dirac-mode number) 2 . On the L 4 lattice, the Dirac-modc 
number is L 4 x N c x 4, which can be reduced to be L 4 x N c , using the Kogut-Susskind technique [2, 4]. Even for the 
projected operators, where the Dirac-mode space is restricted, the matrix is generally still huge. At present, we use 
a small-size lattice in the actual lattice QCD calculation. 

We use SU(3) lattice QCD at = 5.6 (i.e., a ~ 0.25fm) on 6 4 at the quenched level. We show in Fig. 6(a) the 
spectral density p(X) of the QCD Dirac operator p. The chiral property of p leads to p(— X) = p(A). Figure 6(b) is 
the IR-cut Dirac spectral density pm(A) = p(A)6>(|A| — Air) with the IR-cutoff Air = 0.5a _1 ~ 0.4GeV. By removing 
the low- lying Dirac modes, the chiral condensate is largely reduced as {qq)A m /(qq) — 0.02 around the physical region 
ofm- 5MeV. 

Figure 7 shows the inter-quark potential V P (R) after removing low-lying Dirac modes, which is obtained in lattice 
QCD with the IR-cut of pm(A) = p(A)0(|A| - A m ) with the IR-cutoff A m = 0.5a- 1 . Remarkably, no significant 
change is observed on the inter-quark potential besides an irrelevant constant, that is, the confining force is almost 
unchanged, even after removing the low-lying Dirac modes, which are responsible to chiral symmetry breaking. This 
result indicates that one-to-one correspondence does not hold for confinement and chiral symmetry breaking in QCD. 

We also investigate the UV-cut of Dirac modes in lattice QCD, and find that the confining force is almost unchanged 



by the UV-cut. This result seems consistent with the pioneering lattice study of Synatschke-Wipf-Langfeld [2l|: they 
found that the confinement potential is almost reproduced only with low-lying Dirac modes, using the spectral sum of 
the Polyakov loop [22j ■ Furthermore, we examine "intermediate-cut" , where a certain part of Ai < | A„ | < A2 of Dirac 
modes is removed, and obtain almost the same confining force. Then, we conjecture that the "seed" of confinement 
is distributed not only in low-lying Dirac modes but also in a wider region of the Dirac-mode space. 
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FIG. 6: (a) The Dirac spectral density p(A) in lattice QCD at /3=5.6 and 6 4 . The volume V is multiplied, (b) The IR-cut 
Dirac spectral density pm(A) = p(A)0(|A| — Air) with the IR-cutoff Air = 0.5a _1 ~ 0.4GeV. 
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FIG. 7: The circle symbol denotes the lattice QCD result of the inter-quark potential after removing low-lying Dirac modes, 
obtained with the IR-cut of pm(A) = p(A)#(|A| — Air) with the IR-cutoff Air = 0.5a _1 . The square symbol denotes the 
original inter-quark potential. The potential is almost unchanged even after removing the low-lying Dirac modes, apart from 
an irrelevant constant. 



G. Discussions on "Confinement =^ Chiral Symmetry Breaking" in QCD 



Here, we discuss the obtained result of "confinement ^ chiral symmetry breaking" . As for their close relation via 
monopolcs discussed in Sec. 3. 2, the monopole would be essential degrees of freedom for most nonperturbative QCD: 
confinement fl8j . chiral symmetry breaking 15 1, and instantons [2fj| . In fact, removing the monopole would be "too 



fatal" for the nonperturbative properties, so that nonperturbative QCD phenomena are simultaneously lost by their 
cut. 

As shown in Fig. 8, the Dirac zero-mode associated with an instanton is localized around it 0. However, the 
localized objects are hard to contribute to the large-distance phenomenon of confinement, although such low-lying 
Dirac modes contribute to chiral symmetry breaking. 

Recall that instantons contribute to chiral symmetry breaking, but do not directly lead to confinement Q. Then, 
as a thought experiment, if only instantons can be carefully removed from the QCD vacuum, confinement properties 
would be almost unchanged, but the chiral condensate is largely reduced, and accordingly some low-lying Dirac modes 
disappear. Thus, in this case, confinement is almost unchanged, in spite of the large reduction of low-lying Dirac 
modes. 

If their relation is not one-to-one, richer phase structure is expected in QCD. For example, the phase transition 
point can be different between dcconfincmcnt and chiral restoration in the presence of strong electro-magnetic fields, 
because of their nontrivial effect on chiral symmetry [23j . 
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FIG. 8: Around each instanton, the Dirac zero-mode is localized, and such low-lying Dirac modes contribute to chiral symmetry 
breaking. However, the localized objects are hard to contribute to confinement. 



First, we have studied the relevant gluon-momentum region for confinement in lattice QCD, based on the Fourier 
expansion. Remarkably, the string tension a is almost unchanged even after removing the high-momentum gluon 
component above 1.5GeV in the Landau gauge. We then have concluded that confinement originates from the low- 
momentum gluon component below 1.5GeV, which is the upper limit to contribute to a. 

Second, we have developed a manifestly gauge-covariant expansion using the eigen-mode of the QCD Dirac operator 
Tj) = "/^D^. With this method, we have performed a direct investigation of correspondence between confinement and 
chiral symmetry breaking in lattice QCD. As a remarkable fact, the confinement force is almost unchanged even 
after removing the low-lying Dirac modes, which arc responsible to chiral symmetry breaking. This indicates that 
one-to-one correspondence does not hold for between confinement and chiral symmetry breaking in QCD. 
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